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Abstract. The combination of ground-based astromet- 
ric compilation catalogues, such as the FK5 or the GC, 
with the results of the ESA Astrometric Satellite HIP- 
PARCOS produces for many thousands of stars proper 
motions which are significantly more accurate than the 
proper motions derived from the HIPPARCOS observa- 
tions alone. In Paper I (Wielen et al. 1999, A&A 347, 
1046) we have presented a method of combination for 
single stars (SI mode). The present Paper II derives a 
combination method which is appropriate for an ensemble 
of 'apparently single-stars' which contains undetected as- 
trometric binaries. In this case the quasi- instantaneously 
measured HIPPARCOS proper motions and positions are 
affected by 'cosmic errors', caused by the orbital motions 
of the photo-centers of the undetected binaries with re- 
spect to their center-of-mass. In contrast, the ground- 
based data are 'mean values' obtained from a long period 
of observation. We derive a linear 'long-term prediction' 
(LTP mode) for epochs far from the HIPPARCOS epoch 
Th ~ 1991.25, and a linear 'short-term prediction' (STP 
mode) for epochs close to T H . The most accurate predic- 
tion for a position at an arbitrary epoch is provided by 
a smooth, non-linear transition from the STP solution to 
the LTP solution. 

We present an example for the application of our 
method, and we discuss the error budget of our method for 
the FK6 (a combination of the FK5 with HIPPARCOS) 
and for the combination catalogue GC+HIP. For the ba- 
sic fundamental stars, the accuracy of the FK6 proper 
motions in the LTP mode is better than that of the HIP- 
PARCOS proper motions (taking here the cosmic errors 
into account) by a factor of more than 4. 
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1. Introduction 

In Paper I (Wielen et al. 1999b), we have shown that the 
combination of the data of the HIPPARCOS astromet- 
ric satellite (ESA 1997) with ground-based results (such 
as the FK5) is able to provide for many stars individual 
proper motions which are significantly more accurate than 
the HIPPARCOS proper motions themselves. The method 
has been already successfully applied in the construction 
of the FK6, the Sixth Catalogue of Fundamental Stars 
(Part I of the FK6: Wielen et al. 1999d, Part III of the 
FK6: Wielen et al. 2000a). 

The method of combination presented in Paper I is, 
however, strictly valid for single stars only. In the FK6, 
we call this procedure therefore the 'single-star mode'. 

In reality most of the stars are members of binaries or 
of multiple systems. If the duplicity of an individual object 
is already definitely known, either from ground-based in- 
vestigations or from HIPPARCOS observations, then the 
method of combination has to be changed properly in or- 
der to obtain meaningful results. We call such procedures 
'special solutions'. Paper III of this series of papers will 
discuss the special solutions for visual binaries and other 
types of double stars. 

The present paper (Paper II) describes an appropriate 
method of the combination of HIPPARCOS results with 
ground-based observations for 'apparently single stars'. 
Even if we have removed from such a sample of appar- 
ently single stars all the objects with known duplicity, 
there remains beside the truly single-stars a large num- 
ber of hitherto undetected astrometric binaries. The mea- 
sured photo-center of an unresolved astrometric binary 
moves on the sky on a wavy curve, in contrast to the 
linear motion of single stars (Fig. 1). In such a case, an 
'instantaneously' measured proper motion deviates from 
a 'mean' proper motion, averaged over a long interval of 
time. The proper motions provided by HIPPARCOS (ESA 
1997) are essentially such instantaneously measured quan- 
tities, since they are derived from positional measurements 
spread over about three years only. In contrast, the proper 
motions given in the FK5 (Fricke et al. 1988, 1991) are 
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Time t 

Fig. 1. Wavy motion of the photo-center of an astrometric 
binary (solid curve) around the linear motion of its center- 
of-mass (cms; dashed line). An instantaneous position and 
an instantaneous proper motion are indicated (filled dot; 
arrow). The linear prediction based on the instantaneous 
values is shown as a dotted line. The long-term averaged 
'mean' proper motion is equal to the motion of the cms. 

long-term averages over up to about 200 years. We have 
called the difference between the instantaneous proper- 
motion and the mean one the 'cosmic error' of the instan- 
taneous proper motion (Wielen 1995a, b, 1997, Wielen et 
al. 1997). 

In some cases, the individual cosmic error is so large 
that the duplicity of an apparently single star is strongly 
indicated by this fact. We have called such objects 'A^t 
binaries' (Wielen et al. 1999a). In most cases, however, 
the cosmic error in the HIPPARCOS proper motion of a 
hitherto undetected astrometric binary is not individually 
significant, but can be shown to exist only statistically in 
a larger sample of apparently single stars. Nevertheless, 
an appropriate combination method should not neglect 
the statistical consequences of these cosmic errors. Our 
comparison of HIPPARCOS proper motions with ground- 
based results has shown that, at least for brighter stars, 
the average cosmic error in a HIPPARCOS proper mo- 
tion is often larger than the HIPPARCOS measuring er- 
ror, typically by a factor of three (Wielen 1995a, b, Wielen 
et al. 1997, 1998, 1999c). 

Wielen (1997, henceforth called Paper WPSA) has de- 
veloped a coherent scheme of 'statistical astrometry' for 
handling the effects of cosmic errors in high-precision as- 
trometry. In the following sections, we shall apply the con- 
cepts of statistical astrometry to the problem of combin- 
ing the HIPPARCOS results with ground-based measure- 



ments for a sample of apparently single stars. The main 
results will be solutions which we call the 'long-term pre- 
diction mode (LTP)' and the 'short-term prediction mode 
(STP)'. Our results have already been applied for the LTP 
and STP solutions given in Part I and Part III of the FK6 
(Wielen et al. 1999d, 2000a). 

It is clear that in principle the most desirable solu- 
tion for our problem would be to treat each star individ- 
ually and fully correctly. This would mean (1) for truly 
single stars: to use the 'single-star mode (SI)' described 
in Paper I, and (2) for binaries: to apply individual or- 
bital corrections, as e.g. done for Polaris by Wielen et al. 
(2000b). However, since for 'apparently single stars' the 
true nature (single or double) of the individual objects is 
unknown, we have to rely on statistical methods in order 
to handle such a sample of stars properly. The treatment 
of a sample of apparently single stars by our statistical 
procedures gives on average the best astrometric predic- 
tion, and it provides the most realistic error budget for 
such a sample. In this sense, our statistical treatment is 
certainly much more appropriate than to ignore the bi- 
nary nature of a large fraction of a sample of apparently 
single stars altogether. 

In Paper I, we have already pointed out that the older 
observations carry a high weight in the combination of 
ground-based measurements with HIPPARCOS data, and 
that therefore the GC (Boss et al. 1937) should be also 
considered here because of its high number of rather well- 
measured stars. We call the result of the combination 
of the GC with HIPPARCOS the combination catalogue 
GC+HIP. 

2. Basic concepts and equations 

2.1. Concepts 

With respect to the two catalogues which should be com- 
bined into a new one, we follow closely the situation de- 
scribed in Sect. 2 of Paper I. We are using also as far as 
possible the nomenclature of that section. 

We assume that two astrometric compilation cata- 
logues are available, identified by the indices 1 (e.g. for 
the FK5) and 2 (e.g. for HIPPARCOS). For the combined 
catalogue (e.g. the FK6), we use the index C, usually now 
supplemented by an additional subindex which identifies 
the special mode of the solution (e.g. LTP or STP). Each 
of the two basic catalogues (i = 1, 2) provides for the stars 
a position Xi(Ti) at a central epoch Ti and a proper mo- 
tion /ii for two angular coordinates (e.g. right ascension 
a* = a cos 5 and declination S). The mean measuring er- 
rors of Xi(Ti) and /ii are denoted by e x ^ and e^. Usually 
one of the catalogues (e.g. i = 1) has first to be reduced to 
the astrometrical system of the other catalogue (e.g. the 
FK5 to the HIPPARCOS system). In this case, xi and 
\x\ denote the already systematically corrected quantities, 
and their mean errors e Xt i and e^.i include the uncertainty 
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of the systematic corrections. For the determination of 
the systematic differences between two catalogues, we use 
methods developed for the construction of the FK5 (Bien 
et al. 1978). 

With respect to the principles of statistical astrometry, 
we make now the important additional assumptions that 
Catalogue 1 gives 'mean' quantities for x and /i, averaged 
over a long interval of time, while Catalogue 2 provides 
'instantaneously' measured values of x and fx. If we apply 
our scheme to a combination of the FK5 with HIPPAR- 
COS, both assumptions are fulfilled to a high degree of 
approximation: The FK5 is based on ground-based obser- 
vations spread over about two centuries, while the FIIP- 
PARCOS results are obtained from measurements made 
during a short interval of time, about three years only. In 
the terminology of statistical astrometry, our assumptions 
mean that we suppose that the FK5 is free from cosmic 
errors. The cosmic errors in the HIPPARCOS positions 
and proper motions are denoted by c x = (((O)) 1 ^ 2 and 
Cfi = (^(O)) 1 / 2 . The correlation functions ((At), r)(At), 
and ((At) are explained in WPSA (especially Sect. 3). 
Numerical values for c x (p) and c fl (p) as functions of the 
parallax p will be provided in Sect. 5. 

2.2. Basic equations 

In the combination of two astrometric catalogues of which 
at least one contains instantaneously measured data (af- 
fected by cosmic errors), the predicted position x p (t) at an 
arbitrary epoch t should not be anymore a linear function 
of time. According to the principles of statistical astrom- 
etry, the 'best' prediction x p (t) for the true instantaneous 
position x(t) is given by the non-linear expression (WPSA, 
Eq. (82), with a slight change in nomenclature): 

x p (t) = ari(Ti)+^i(t-Ti) 

+ l{t) (x 2 (T 2 ) - a;i(Ti) - m (T 2 - T x )) 

+ (3(t) (MT 2 ) - Mi) (* - T 2 ) . (1) 

The quantities j(t) and (3(t) are functions of time t to be 
determined. For computational reasons, we often use 

B(t)=0m-T 2 ) (2) 

instead of (3(t). 

The functions 7(f) and [3(t) are determined from the 
condition that the mean error s x . p (t) of the predicted po- 
sitions x p (t) should be a minimum for every value of t, 
averaged over the ensemble: 

< P (*) = < 0*p(*) - x (t)) 2 > = ™n. ■ (3) 

The operator < q > means, as in WPSA (1997), the av- 
erage of the quantity q over the ensemble of similar stars. 
For an individual star, < q > can be interpreted as the 
'expectation value' of the quantity q for this star. In this 
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sense, we consider the ensemble averages as statistical pre- 
dictions for the mean behaviour of a typical individual 
member of the ensemble. 

Inserting Eq. (1) into Eq. (3), and using the scheme of 
statistical astrometry, we obtain for e xp (t): 

<p(*) =< ((*p(*)-*m(t)) - (x(t) - Xm(t))) 2 > 

= < ((a:i(T 1 )- a5m (T 1 ))(l-7(t)) 

+ (tn - ihn){t - Ti - 7 (t)(T 2 - Ti) - /?(t)(t - T 2 )) 
+ 7(t)((z 2 (T 2 ) - x(T 2 )) + (x(T 2 ) - x m (T 2 ))) 
+ /3(t)(t - T 2 )((ii 2 (T 2 ) - n(T 2 )) + (m(T 2 ) - Mm)) 

- (x(t) - x m (t)fj > 
= (l + (7W) 2 K(0)-2 7 (iK(t-T 2 ) 

+ (f3(t)(t - T 2 ))\(0) -2p(t)(t - T 2 )((t - T 2 ) 

+ <i(l-7(i)) 2 +< 2 (7(*)) 2 

+ -Ti- 7 (t)(T 2 - Ti) - p(t)(t - T 2 )) 2 

+ el. 2 (P(t)(t-T 2 )f . (4) 

The quantities x m (t) and fi m — x m are the true mean 
position and the true mean proper motion of the star (in 
the sense of statistical astrometry) . x m is a linear function 
of t, e.g. 

x m (t)=x m (T 1 ) + t i m (t-T 1 ) . (5) 

In deriving Eq. (4) we have used the following relations: 
According to our assumption that Catalogue 1 provides 
mean quantities, xi(Ti) and \x\ differ from x m (Ti) and 
[i m by their measuring errors only: 

<( a;i (r 1 )-x m (r 1 )) 2 > = e 2 4 , (6) 

< (Pi - Mm)' > = 4l • (7) 

The quantities x(t) and /i(t) are the true instantaneous 
position and proper motion of the star at epoch t. Since 
Catalogue 2 is assumed to provide instantaneous quanti- 
ties, x 2 (T 2 ) and fi 2 (T 2 ) differ from x(T 2 ) and fJb(T 2 ) also 
by their measuring errors only: 

< (x 2 (T 2 ) - x(T 2 )f > = el a , (8) 

< (P2(T 2 ) - KT 2 )f > = < 2 • (9) 

The correlation functions £(At), n(At), and ((At) (see 
WPSA) occur because of: 

< (x(T 2 ) - x m (T 2 )) 2 > = < (x(t) - x m (t)) 2 > = ((0) ,(10) 

< (x(T 2 ) - x m (T 2 ))(x(t) - x m (t)) > = ((t- T 2 ) , (11) 

< (fi(T 2 ) - /i m ) 2 > = ,7(0) , (12) 

< (x(t) - x m (t))(fi(T 2 ) -n m )>= ((t - T 2 ) . (13) 
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We have further to remember 



with the auxiliary quantity N: 



< (x(T 2 ) - x m (T 2 ))(p(T 2 ) - Mm ) > = C(0) = . (14) 

The products of the other terms occuring in the second 
line of Eq. (4) all vanish, because the measuring errors 
are not correlated with the cosmic errors, for example 



N 



< (a;i(Ti) - x m (T!))(x(t) - x m (t)) > = , 



(15) 



and because the measuring errors of Xi(Ti) and Hi{Ti) are 
not correlated at the central epoch Tj (by definition of the 
central epoch), for example: 



< (^(Ti) - a; m (Ti))(/ii - n m ) > = 



(16) 



and finally because the measuring errors of x\ and \i\ are 
not correlated with those of x 2 and fi 2 , for example 



< (a:i(Ti)-a; m (Ti))0i 2 (T 2 )-/i(T2)) > = . 



(17) 



The conditions for the minimum of e^, p (t) with respect 
to 7 and /3 are 



94 P 

and 

del v 

— ^ = . 

dp 



(18) 



(19) 



Carrying out these procedures and using the final part of 
Eq. (4) for e^ p (i), we obtain the following two equations 
of condition for the functions j(t) and (3{t) which give the 
'best' values of s XtP (t): 

7(*)(£(0) + 4,1 + 4,2 + £ l,i( T 2 - T i) 2 ) 
+ /3(t)(t-T 2 ) £ 2 il (T 2 -T 1 ) 

= e(t-T 2 )+4 1+e 2 jl ( t _T 1 )(T 2 -T 1 ) , (20) 

7 (i)4, 1 (T 2 -T 1 ) 
+ /J(f)(t-T 2 )(r ? (0)+4 il +£^ 2 ) 

= C(t-T 2 )+^ il (t-Ti). (21) 

The Eqs. (20) and (21) correspond to the Eqs. (86) and 
(87) in WPSA. Solving the Eqs. (20) and (21) for the 
unknowns j(t) and (3(t) we obtain: 

7(*) - ((£(* - r 2 ) + <i)(»?(o) + 4,! + 4, 2 ) 

+ 4,i( r 2 - TiK^O) + 4 2 )(i - Tr) - C(t - T 2 )))/J\f(22) 

fl(t) = /3(t)(t-T 2 ) 

= (C(* - T 2 )(£(0) + 4,1 + 4,2 + - ?i) 2 ) 

+ 4i( i - T i)(^°) + 4,i + 4, 2 ) 

- 4i(r 2 - 7i)(£(i - T 2 ) + 4,i))/^V , (23) 



2 , 2 

x,l ^ G :r 



2)(^(o) + 4i + 4, 2 ) 



ri) 2 (^(o) + 4 2 ). 



(24) 



For deriving the functions j(t) and /3(t), we have up 
to now implicitely assumed that all the stars in the en- 
semble have the same measuring errors. This is, of course, 
not strictly true in reality. However, the ensemble averages 
are actually neccessary for handling the cosmic errors only, 
but not for treating the measuring errors. Hence we shall 
use all the equations derived above by inserting the indi- 
vidual measuring errors if we treat individual stars of the 
ensemble. Formally we may imagine to handle subsamples 
of stars in which the stars have the overall behaviour with 
respect to the cosmic errors, but in which the common 
measuring errors are equal to those of the individual star 
under consideration. It is a more severe problem that we 
do not make use of our knowledge of how large the indi- 
vidual cosmic errors in x and fi are for a given individual 
object (e.g. [i 2 — m). As discussed in WPSA, this would 
require 'conditioned' correlation functions. Since this in- 
formation is presently not available, we are treating here 
the consequences of the cosmic errors on the level of en- 
semble averages only. 

Inserting these results for 7(f) and (3(t) from Eqs. (22) 
and (23) into Eqs. (1) and (4), we derive the prediction 
x p (t) for the instantaneous position x(t) of the star and 
the mean error £ x ,p{t) of this prediction. The prediction 
x p (t) is a non-linear function oft, because the correlation 
functions £,{t — T 2 ) and ((t—T 2 ), which occur in the formu- 
lae for 7(t) and 0(t), are non- linear functions. A typical 
run of x p (t) is shown in Fig. 2. 

In order to illustrate the properties of our prediction 
Xp(t), we consider in the following Subsections 2.3 and 2.4 
two limiting cases in which we either neglect the measuring 
errors or the cosmic errors. 



2.3. Measuring errors neglected 

If we neglect all the measuring errors and set £ x ,i(Ti) = 
£x 2(^2) = £ p 1 = e M 2 =0, then we obtain for -f(t) and 
0(t) 

7nmc(t) = £(t - T 2 )/£(0) , (25) 
B nmc (t) = /3 nmo (t)(t - T 2 ) = ((t - T 2 )/r)(0) . (26) 

These results were already derived and discussed in WPSA 
(Sect. 4.2.4 and Fig. 10). The corresponding prediction 
x Pl nme(t), shown in Fig. 2, passes through the point x 2 (T 2 ) 
with the slope fi 2 (T 2 ). Therefore, the measured instanta- 
neous position and proper motion at epoch T 2 are ex- 
actly reproduced by x p , nmc (t). For t — > ±00, the predic- 
tion x p ^ nmc (t) approaches asymptotically the mean posi- 
tion x m (t) = xi(t) = x 1 (T 1 ) +fj, 1 (t-T 1 ) of the star. The 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



5 



uncertainty £ x ,p,nme(i) of the prediction is given by 
(£(0)) 2 -(£(t-T 2 )) 2 {C{t-T 2 )f 



'x,p,nme 



it) 



7,(0) 



(27) 



We see that £ l p „ mc is zero at the epoch T 2 and approaches 
c x = (^(O)) 1 ^ 2 for i — > ±oo, since in this limit a; p (i) is 
equal to x m (t), and ir m (i) differs from the instantaneous 
position x(t) by the cosmic error c x on average. 

2.4- Cosmic errors neglected 

If we neglect the cosmic errors and set £ = n = ( = 0, 
then we obtain 



7. 



>(*) — ( e s,l( £ f*,l + £ f*>2) 

+ 4, 1 4 2 (T 2 -r 1 )(i-T 1 ))/Ar nce , (28) 

Snce(i) = /3nce(i)(i - T 2 ) = ( - 4,l^,l( T 2 ~ T l) 

+ (4i + 4 2 )<i(^-r 1 ))/iV ncc , (29) 
with 



( £ L+ £ L)( £ u,l+ £ L) + 



>(T 2 -T 1 ) 2 .(30) 



Equation (29) illuminates one of the advantages of in- 
troducing B(t) — [3(t)(t — T 2 ) as a substitute for (3(t). 
While B nce (T 2 ), and in general B(T 2 ) for non-zero measur- 
ing errors, remains finite, the quantity (3 ncc , and in general 
j3, tends towards infinity for t — > T 2 . Only in some degen- 
erated cases, such as £„ i = or e Xi2 = or e x \ — > oo, the 
quantity /3 n ce 5 and in general /3, remains finite for t —> Ti- 
ll we insert 7 nC e(^) and B nce (t) into Eq. (1), the corre- 
sponding prediction x p ^ ncc (t) for the position of the star 
at an epoch t is now a strictly linear function of t, since 
7ncc(t) and B ncc (t) are linear in t. 

In order to facilitate the understanding of the be- 
haviour of the prediction x PtIlce (t), we rewrite a; Pincc by 
using the auxiliary quantities Tc,ncc, £c,ncc, Mc,ncc- If we 
insert Eqs. (28) and (29) into Eq. (1), we obtain after some 
lengthy algebra: 



(Tc,nce) + MC.nco {t - Tq 

with the auxiliary quantities 



T. 



C.nce 



^C.nce (Tb,nce) = 



^C,ncc 

Ho 



w Xi iT x + w x , 2 T 2 
w x ,iXi{T\) + w Xt2 x 2 (T 2 ) 

W x .l + W x . 2 

(x 2 (T 2 )-x 1 (T 1 ))/(T 2 -T 1 ) 



(31) 

(32) 
(33) 

(34) 
(35) 



The weights w are given by 
1 

Wx,l = ~2— . 

£ x,l 



W x2 



«V,1 



-x.2 



-fi,2 



1 (T 2 - Tr) 2 



+ £ x,2 



(36) 
(37) 
(38) 
(39) 
(40) 



Inserting 7 nC c(i) and B nce (t) into Eq. (4), we obtain the 
mean error e x . p ^ ncc (t) of x Pjnco . Using the form of x Pincc 
as given in Eq. (31) and the auxiliary quantities with the 
index C, we find 

£ x,p,nce(*) — £ lancc(^C,ncc) + £ ^,C,nce (* ™ ^C,nco) 2 , (41) 

with 

1 



e x,C,nce(^C,nce) — 



Wx.l + W Xi2 



1 



-"^,C,nce 



Wn,i + w^ 2 + w^o 



(42) 



(43) 



A comparison of the Eqs. (31)-(42) with the analytic 
version of the single-star-mode solution of Paper I (Sect. 
2, especially Eqs. (19), (23), (30)-(32)) shows that the pre- 
diction 

^p,ncc(^) &nd its mean error £aj,p,nce 

(t) are strictly 

identical with the single-star- mode solution xsi{t) and its 
mean error e Xl si(i)- This result is very pleasing, since it 
proves the internal consistency of our scheme: In the limit 
of vanishing cosmic errors, the prediction x p (t) according 
to Eq. (1) is asymptotically approaching the single-star 
mode solution xsi(t). This result is not apriori self-evident, 
since our definitions for the 'best' solution for predicting 
x(t) differ at least formally in Paper I and in this Paper II 
(i.e., Eqs. (5) and (35) of Paper I versus Eq. (3) of Paper 
II). 

2.5. Motivation for introducing the long-term and 
short-term prediction 

The general solution of our combination problem is given 
in Sect. 2.2. for all epochs t. The solution x p (t) is a non- 
linear function of t, and requires the knowledge of the 
correlation functions £(At) and ((At) as functions of the 
epoch difference At. At present, we do not have a well- 
established knowledge about the run of £(At) and C(Ai). 
Only the cosmic errors c p = (r/(0)) 1//2 and c x = (C(O)) 1 / 2 
can be empirically determined from the comparison of the 
FK5 with HIPPARCOS, assuming that the FK5 is giving 
'mean' quantities. 
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Even if we would know the run of £(At) and C(Ai) as 
a function of At, the non-linearity of x p (t) would demand 
a table of x p (t) for a sequence of epochs, e.g. for each 
year, if the user should not have the burden to do the full 
calculation himself by running a program. 

We propose the following solution: The general solu- 
tion for x p (t) allows rather easily to obtain two limiting 
solutions for At = t — T 2 — > ± oo and for At — > 0. We 
call the solution for At — > ± oo the 'long-term prediction 
(LTP)', and the solution for At — > the 'short-term pre- 
diction (STP)' around the epoch T 2 . 

Both the LTP and STP solutions are linear in t. They 
can be therefore given in the usual astrometric style, i.e. 
as a position at a central epoch and a proper motion, 
together with their mean errors. The details on the LTP 
and STP solutions xltp(^) and xstpW are given in the 
Sects. 3 and 4. 

The general solution x p (t) is a smooth transition from 
short-term prediction xstpW (for epochs around T 2 ) to 
the long-term prediction xltp(0 for epochs t far away 
from T 2 . In Sect. 5, we will discuss a convenient (but only 
approximately valid) method to carry out this transition, 
if we know the run of C(At). This gives at least a rough 
indication on the process of transition as a function of the 
epoch difference At, even if C(At) is not well-established. 
If, in the future, C(At) should be better determined, then 
our method would allow rather conveniently the (approx- 
imate) determination of x p (t) also for epochs inbetween 
the validity ranges of xstp(£) and Xltp(£)- 



3. Long-term prediction (LTP) 

We consider in this section the limiting case of the gen- 
eral solution x p (t) for \t — T 2 \ — > ± oo. This 'long-term 
prediction' xltp(^) is valid for epochs not too close to the 
epoch T 2 of the instantaneous Catalogue 2 (i.e. here the 
HIPPARCOS Catalogue with T 2 - 1991.25). 

We assume that the epoch difference At = t — T 2 is so 
large that the correlation functions £(At) and C(At) both 
vanish. Setting £(t — T 2 ) = and ((t — T 2 ) = 0, we obtain 
from the general Eqs. (22) and (23) for the LTP solution 

7ltpW - (4,1 (4,1+ K,2 + »7(0)]) 

+4,1 [4,2 + V(0)] (T 2 - Tr) (t - Tr) ) /N , (44) 
B LT p(t) = Pltp (t)(t-T 2 ) 

= ( - 4i4,i ( T 2 - to 

+ (4,1 + [4,2 + £(0)] ) 4,i (t - Tr) ) /N , (45) 

where N is still given by Eq. (24). If we compare the Eqs. 
(43), (44), and (24) for 7ltp, -B LT p, and N = N LT p with 
the corresponding Eqs. (28)-(30) for 7 nce , B ncc , and N nce , 



we find that they are identical if we replace in the equa- 
tions for the nee solution the quantity £ 2 X 2 by 



42, LTP ~ £ ^,2 + ^(0) — £ x,2 + C . 

and £^ 2 by 



,2 

x ' 



e M ,2,LTP — £ n,2 



-2 

-M,2 



4- 



(46) 



(47) 



This is very plausible, since the instantaneously measured 
quantities x 2 (T 2 ) and fi 2 (T 2 ) are affected by the cosmic 
errors c x and c M . If we add these cosmic errors to the 
corresponding measuring errors £ x , 2 and £^ 2 , we obtain 
'apparent' measuring errors £ Xi2 ,ltp and £ Mj2 .ltp- Since 
the cosmic errors are not correlated with the measuring 
errors, the summation has to be done quadratically. 

Using this finding we obtain for the long-term predic- 
tion 



xltpW = x LT p(T LT p) +MLTP (*-T LT p) , (48) 



with 



Tltp = 



xltp(?ltp) — 



Mltp 



W x ,l T\ + Wa,2,LTP T 2 

Wx,l + W X ,2,LTP 
U>z,l Xi(T'i) +W I[2 ,LTP x 2 (T 2 ) 

Wx,l + Wa;,2,LTP 
«V,1 Ml + W Mi2 ,LTP M2 + WVALTP Mo 



Wf+.i + W M , 2: LTP + «V,o,LTP 

The weights w are given by Eqs. (36) and (38), and by 
1 



(49) 
(50) 
(51) 



11^,2, LTP 



£ x,2 + C x 



1 



«V,2,LTP - -J- 



Wp,o,LTP — — 



(T 2 - Ti) 



-x.l 



+ £ x,2 + c x 



(52) 
(53) 
(54) 



The long-term prediction Xltp (t) has in fact two con- 
ceptionally different properties: (1) As described above, 
it is the limit of the general prediction x p (t) for the true 
instantaneous position x(t) for t — > ± oo. (2) On the other 
hand, xltp(^) is for all epochs t the best prediction for 
the mean position x m (t) of the object. This means espe- 
cially that mltp is the best estimate of the center-of-mass 
velocity of the object. 

The two different concepts produce two different error 
estimates £ x ,ltp(£) for xltp(^)- If we consider xltp as 
the prediction for the mean position x m (t) then the mean 
error £ x ,LTP,m(*) is given by: 

e x,LTP,m(*) = 4LTP,m(7LTp) + 4,LTP,m (* _ ?LTp) 2 ,(55) 

with 



4,LTP,m(-^LTp) — 



1 



w 



x.l 



w 



(56) 



£c,2,LTP 
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-jli.LTP.m 



(57) 



>,2,LTP 



>,0,LTP 



If we consider xltp as the prediction for the instantaneous 
position x(t) for large values of |At| = \t — T 2 \, then the 
uncertainty £ x ,LTP,in S (*) of x LTP (i) is given by 



£ 2,LTP,ins(*) — £ x,LTP,mW + c x ■ 



(58) 



Similarly, the mean error e^.LTP.ins of the predicted in- 
stantaneous proper motion for large At is given by 



2 2 , 2 

£ M,LTP,ms — £^,LTP,m + C fj, 



(59) 



The equations above describe what we have called the 
'analytic' approach in Paper I. For the 'numerical' ap- 
proach, we can take over for the LTP solution the for- 
mulae given in Paper I for the single-star mode (SI) with 
the following changes: (1) We should not redetermine the 
parallax p of the star. The determination of p by HIP- 
PARCOS requires instantaneous values of x and p,, not 
the 'mean' LTP values. Formally, we set all correlations 
between p and the other quantities (position, proper mo- 
tion) equal to zero. (2) In the variance-covariance matrix 
D we replace in the diagonal line the HIPPARCOS values 
of Alii = 4*,h(?h) by 



-Dh,ii,ltp - sL.hPh) + c l 



(60) 



and e| H by e| H + c 2 , s 2 ^^ by e\ a , + c 2 , and e^ H 
by e 2 ^ s H + c 2 , respectively. No changes are made in the 
non-diagonal elements (except for decoupling the parallax 
as described above), i.e. all the covariances remain in the 
LTP as they were in the SI mode, e.g. 

£>h,12,ltp = -Dh,i2 = e a *,n{Tn)£s,ii(Tii)pa6,ii(Tu) ■ (61) 



4. Short-term prediction (STP) 

In this section, we consider the other limiting case of 
the general solution x p (t), namely the limit for At = 
t — T 2 — > 0. This 'short-term prediction' xstpW is valid 
for epochs close to T 2 (in the case of using HIPPARCOS: 
T 2 ~ 1991.25). 

For the case At — > 0, we use for the correlation function 
£(Af) and ((At) Taylor series in At, and keep only terms 
which are linear in At. From the Eqs. (54) and (55) of 
WPSA, we obtain for small values of At = t — T 2 



£(At) 
C(At) 



n(0)At 



(62) 
(63) 



In Eq. (63), we have made use of the differential rela- 
tion (28) of WPSA. Although HIPPARCOS values are 
already averaged over about 3 years of observation, the 
use of Eq. (63) is justified according to numerical investi- 
gations carried out by M. Biermann (1996, unpublished, 
see also WPSA, Sect. 3.4). 



Inserting Eqs. (62) and (63) into the Eqs. (22) and 
(23), we find for the short-term prediction 

7STp(t) = ((£(0) + 4,i)(»?(0) + 4i + 4,2) 
+ 4i( T 2-ri) 2 (ry(0) + 4 2 ) 

+ 4, 1 4, 2 (T 2 -T 1 )(t-T 2 ))/AT, (64) 
Ssxp(i) = Pstp (t)(t-T 2 ) 

+ (v(0) (£(0) + 4,i + 4.2 + elAT 2 - T1) 2 ) 

+ 4i(C(0)+4 1+ e 2 , 2 ))(t-r 2 ))/iV, (65) 

where TV is still given by Eq. (24). Inserting 7stp(*) and 
BsTp(t) into Eq. (1), we obtain 



xsTp(t) = x ST p(T 2 ) + ustp (t-T 2 ) , 



(66) 



which is a linear function of the epoch t. We do not use 
here a 'central' epoch Tgxp, since our Taylor series for 
XsTp(t) is centered around t = T 2 . The proper motion 
Mstp can be derived as the sum of the coefficients in front 
of t in the Eqs. (1), (63), and (64): 

Mstp = Mi 

+ ((z 2 (T 2 ) - ^(Tr) - Ml (T 2 - 7i)) sl A el 2 (T 2 - T^/N 
+ (M2 - Mi) (v(0)(m + <i + < 2 + el x (T 2 - T,) 2 ) 

+ 4i(£(0)+eL+4 2 ))/AM67) 
After some algebra, mstp can be rewritten as 

«V,10,STP MlO + WpL,2 M2 



Mstp 



"V.io.stp + w^ 2 



(68) 



The 'combined' mean proper motion /i 10 is derived from 
Mi and /z by 



M10 



Mi + "Valtp Mo 

«V,1 + W Mi0 ,LTP 



(69) 



MOj and w^^ are given by the corresponding equa- 

tions in Sect. 2.4, and u^o, ltp by Eq. (53). The weight 
w m ,io,stp of jUio in the STP solution is given by 



1 



«V,10,STP = 



with 



(I/uV.io.ltp) + e^ 10)LT p + c; 



W M ,10,LTP = + UV >0 ,LTP = — 



, (70) 



(71) 



-jMO.ltp 



The meaning of Eqs. (68)-(71) is the following: mstp is the 
weighted mean of /i 2 and the combined mean proper mo- 
tion ni , where the cosmic error c M in /Lt 10 has to be taken 
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into account in w m .io,stp- The combined mean proper mo- 
tion /iiio is itself a weighted mean of [i\ and /j,q, where for 
Ho the cosmic error c x in X2(T 2 ) has to be included into 

U^ALTP- 

The exact expression for xsTp(Tb) is given by 
xstp(T 2 ) = x 2 (T 2 ) 

- 4,2((^(0) + 4i + sl 2 )(x 2 (T 2 ) - nCTi) - Mi(T 2 - Ti)) 

-^(Ta-TO^Ta)-/*!))/^. (72) 

For all practical purposes we can, however, neglect the last 
term in Eq. (72), which is proportional to e 2 x 2 . The mea- 
suring accuracy e Xj2 for the HIPPARCOS positions x 2 (T 2 ) 
at T 2 ~ 1991.25 is so high, relative to the other measuring 
errors and cosmic errors, that this term is usually of the 
order of 0.01 mas only. We use instead of Eq. (72) the very 
good approximation 

xstv{T 2 ) = x 2 {T 2 ) . (73) 
The mean errors of xstp(T 2 ) and of /^stp are given by 
4sTp(r 2 ) - 4, 2 (T 2 ) , (74) 



In presenting the results in printed form, we use again cen- 
tral epochs Tstp (different for a* and 8), at which xstp 
and /xstp are uncorrelated. However, Tstp differs usually 
only very slightly from the individual central epochs of the 
basic HIPPARCOS data. 

5. Transition from the short-term prediction to 
the long-term prediction 



1 



-M,STP 



W M ,10,STP + 10^2 



(75) 



The full uncertainty £ Xj sTp(i) of £stp(0 is given by 

£ x,STp(*) = £ x,STp( 7 2) + e^,STP (t - T 2 ) 2 



(76) 



In Eq. (76), we have neglected the very small correlation 
between Xstp{T 2 ) — x 2 (T 2 ) and ustp which is caused by 
the use of £2(^2) in deriving no which in turn enters into 
Mstp- The last term in Eq. (76) is the statistical uncer- 
tainty of a prediction based on instantaneous data at T 2 
for small epoch differences At (see Eq. (57) of WPSA). 
£,q IV ^ is the fourth derivative of the correlation function 
£(At) with respect to At at At = 0. 

In the numerical approach for deriving the short-term 
prediction, we modify the procedure for the single-star 
mode presented in Paper I more strongly than for the 
LTP solution. As 'observations' b we use first bj^Tn), 
as in Paper I, the parallax pn inclusive. The correspond- 
ing part Dh of the variance-covariancc matrix D remains 
also unchanged. The second part of b, which we now call 
b m is given by the a* and 5 components of the combined 
mean proper motion /zio. The quantity /j,w is obtained in 
a preparatory step from Eq. (69), and its mean error from 
Eq. (70) as (w Mj io,stp) -1 ^ 2 - Each component of /iio, i.e. of 
b m , is considered not to be correlated with any other com- 
ponent of b. The numerical approach for STP produces 
values for x§tp and /^stp (i n cx* and 6), a new parallax 
PSTP: and the corresponding variance-covariance matrix. 



5.1. Transition in position 

As discussed in Sect. 2.4 and illustrated in Fig. 2, the 
general solution x p (t) is a smooth transition from the 
short-term prediction xstp(^) for epochs t close to T 2 to 
the long-term prediction xltp(*) for t — > ±00. We are 
now asking for the 'transition function' /3t ians (t) which de- 
scribes this transition in x: 

X p (t) = (1 - Arans(t)) X LTP (t) + AransW X STP (t) . (77) 

Formally we can always solve Eq. (77) for /3t iana (t), using 
our former results for x p (t), xltp(£), and xstpW- How- 
ever, the resulting transition function Pt ians (t) is then very 
complicated and depends unfortunately not only on the 
correlation functions, but also explicitcly on the measured 
values of x\, x 2 , Mi> M2 ; and on their measuring errors. Such 
a result is not very suitable for a practical application. 

There does exist, however, an approximate treatment 
for the transition function /3t ians (t) which gives a very sim- 
ple and easily applicable result, and which nevertheless 
describes the transition quantitatively rather accurately. 
The basic idea is the observation that in real applica- 
tions the mean position x\(Ti) enters into the final re- 
sult x p (t) mainly through the proper motion fi . This is 
caused by the small error of the HIPPARCOS position 
x 2 (T 2 ) with respect to the error of the mean (FK5 or GC) 
position xi(Ti). Only in cases of a large cosmic error c x {p) 
in x 2 (T 2 ), our approximation becomes less accurate. We 
therefore consider the transition function /3trans(t) for the 
limiting case in which e x ,\ tends towards infinity while 
e^fl remains finite (equal to its actual value). The latter 
can be enforced by setting T 2 — T\ = £a;,i/ e ^,o- This means 
that we let go Ti — > — oo and e Xt i — > + oo in such a way 
that e^fi remains constant. 

If we use this special case as an approximation, we de- 
rive after some lengthy algebra the following rather simple 
expression for the transition function: 



/^trans(^) 



V(0) (t - T 2 ) ■ 



(78) 



A similar function has already been derived as Eq. (70) in 
Sect. 4.2.4 of WPSA. The function /3 t rans has the welcome 
property that it does neither depend on the measured val- 
ues of x and /i nor on the mean errors of a given star. 
These values are fully absorbed in the individual solutions 
xsTp{t) and XLTp(i)- Hence Arans(*) is the same function 
for all the objects. We should remark here that the ap- 
parently more complicated form of the transition function 
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tt — i — r 

/ STP 




t - T 3 [years] 



Fig. 2. Schematic illustration of the non-linear transition 
of the 'best' prediction x p (t) (solid curve) from the lin- 
ear short-term prediction xsTp(i) (STP mode, valid for 
epochs t ~ T 2 ; short-dashed line) to the linear long- 
term prediction xltp(^) (LTP mode, valid for epochs 
\t — T2I — > 00; long-dashed line). 

Table 1. Transition functions /?trans(i) and v trans (t) for 
a simple example 



\At\ = \t-T 2 \ 
[years] 


Arans(At) 


f trana (A£) 





1.000 


1.000 


1 


0.978 


0.934 


2 


0.915 


0.752 


3 


0.819 


0.491 


4 


0.700 


0.202 


5 


0.573 


- 0.065 


6 


0.449 


- 0.270 


7 


0.336 


- 0.397 


8 


0.241 


- 0.445 


9 


0.165 


- 0.430 


10 


0.108 


- 0.373 


11 


0.068 


- 0.297 


12 


0.041 


- 0.219 


13 


0.023 


- 0.152 


14 


0.013 


- 0.099 


15 


0.007 


- 0.060 


20 


0.0001 


- 0.002 



according to Eq. (75) of WPSA is caused by the fact that 
we have used in WPSA a slightly different definition of 
the long and short-term prediction, namely X2(t) instead 
of xstpW and x\(t) instead of xltp(*)- 

The function f3 tTa , ns (t) has the desired properties in the 
limits At — > or 00. For t = T 2 , we have /?trans(^2) = 1) 



because of C(At) ~ »7(0)At for small values of At. For 
t — > 00, ((t — T 2 ) vanishes and hence /3trans(°o) = 0. An 
example for the full run of (3 tr ans(t — T2), from 1 to 0, is 
illustrated by the full curve shown in Fig. 8 of WPSA. Ta- 
ble 1 gives /3trans(i — T 2 ) for a few values of \t — T 2 \, using 
the example given in Sect. 3.6 of WPSA for the correlation 
function C(Ai). At \t — T 2 \ ~ 5-6 years, the general solu- 
tion x p (t) is about half-way between xstp (t) and ccltp (t) ■ 
From both the Fig. 8 of WPSA and Table 1, we get an 
indication for the range of applicability of the short and 
long-term prediction. The short-term prediction xstpW 
has a rather limited range of applicability, namely a few 
years around T 2 only. The long-term prediction xltp(^) 
is a good approximation for the general solution x p (t) for 
epoch differences \t — T 2 \ which are larger than about 10 
years. Hence the transition from STP to LTP is rather 
rapid, at least in our example. 

We should remark here that (by chance) the transition 
function /?trans(i) given by Eq. (78) is even strictly valid 
(exact) for the example used for the correlation functions 
in Sect. 3.6 of WPSA and adopted in Table 1. For more 
general runs of the correlation functions, /3trans can be- 
come quite large at epochs around the 'crossing time' at 
which xltp = a;sTP, if ^ltp and xstp are significantly 
different from x p at that time. 

5.2. Transition in proper motion 

The most accurate prediction fj, p (t) for the instantaneous 
proper motion /j,(t) at an arbitrary epoch t is formally 
given by 

Mp(*) = Xp(t) = (1 - "trans. (*)) A^LTP + "trans (*) MSTP ■ (79) 

The transition function ^trans(^) for the proper motion /z 
is, similar to /3trans, in general a rather complicated func- 
tion, which depends on the measured values of /xi, fi 2 , and 
/j,q. If we adopt the same approximation as used in Sect. 
5.1 for obtaining Eq. (78) for /3trans> namely £ X; i — > oo 
while keeping e^q constant, we derive for ftrans(t) also a 
very simple expression: 

^trans(i) = Arans(t)+Arans(*)(*-7sTp) = ~ "Jf^ -CSQ) 

r?(0) 

This transition function ^trans(i) for n(t) has the same 
proper limiting values as (3trans(t). For t — > oo, we have 
ftrans(co) = and hence ^ p (oo) = /xltp, and for t = Tstp 
we obtain ^ tr ans(7sTp) = 1 and ^(Tstp) = Mstp • In 
Table 1 we list f tra ns(At) for a few values of |At| = \t — 
Tstp |, using Eq. (80) and the simple example for rj(At) 
given in Sect. 3.6 of WPSA. 

The transition functions /3trans(£) and v trans (*) accord- 
ing to Eqs. (78) and (80) have another nice property: If 
we use the example given in Sect. 3.6 of WPSA, then the 
function /3 t rans depends only on t — T 2 , but not on the 
individual cosmic error c ll (p) — (^(O,/))) 1 / 2 , since 77(0) oc- 
curs also as a factor in ((t — T 2 ) and cancels out in Eq. 
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(78). Similarly, the transition function ft ra ns(i) is a 'scaled' 
function and depends on t — T 2 only, but not on r)(0,p). 

6. Cosmic errors 

For deriving the solutions in the LTP and STP mode, 
we need to know the cosmic errors c M = (77(C))) 1 / 2 and 
c x = (C(O)) 1 / 2 . The cosmic errors depend strongly on the 
distance r of the star from the Sun, or equivalcntly on 
the stellar parallax p. For 1202 'apparently single stars' 
from the basic FK5, we have obtained the functions c M (p) 
and c x (p) empirically by using groups of stars in various 
distance intervals. The data can be represented by the 
following fit functions: 

c» = fo(0,p))V» = ( (c| ^ 2)1/2 ) 1/2 , (81) 

c x (p) = (C(0,P)) 1/2 = C 3 c„(p) , (82) 
with 

Ci = 9.30 (mas/year) 2 , (83) 

C 2 = 22.14 mas , (84) 

C 3 = 5.93 years . (85) 

These versions of c M (p) and c x (p) have been used for the 
FK6 (Wiclen ct al. 1999d, 2000a). A table for c M and c x 
as a function of p or r, based on Eqs. (81)-(85), is given 
on page 12 of Wielen et al. (1999d). 

In the future we hope to determine also the depen- 
dence of the cosmic errors on the absolute magnitude (or 
mass) of the stars for a given parallax. A comparison of 
the results for the cosmic errors based on the FK5 stars 
with those derived from the (on average fainter) GC stars 
seems to indicate a rather weak dependence on the bright- 
ness of the stars (Wielen et al. 1998). 

7. An example: a Ari 

In order to illustrate our combination method for all the 
three modes (SI, LTP, STP), we give in Table 2 the results 
for one individual star. As in Paper I (Sect. 5), we give the 
positions x{t) and the proper motions /i(t) always relative 
to the HIPPARCOS solution as Ax(t) = x(t) - x H (t) and 
A^(t) = ^(t) — Hn{t), in order to save printing space and 
to make the comparison of the results easier. The ground- 
based data are always reduced to the HIPPARCOS sys- 
tem. The results of the single-star mode presented here in 
Paper II differ somewhat from those of Paper I, because we 
now adopt slightly improved systematic differences FK5- 
HIP and GC-HIP. For a valid comparison of the results of 
the three different modes it is necessary to use exactly the 
same basic input data. 

Table 2 shows that the short-term prediction is usu- 
ally quite close to the HIPPARCOS solution. On the other 
hand, the long-term prediction differs most strongly from 



the HIPPARCOS solution, since the HIPPARCOS values 
are entering with a lower weight into the LTP mode than 
in the SI mode, because of the cosmic errors in the HIP- 
PARCOS data. The central epoch Tltp is the only one 
which is usually significantly earlier than Th ~ 1991.25, 
and the mean error of the central position £ltp(7ltp) is 
typically only slightly smaller than the cosmic error c x (p). 



8. Error budget 

In Table 3 we present the error budget of proper motions 
in the three different modes (SI, LTP, STP) for two sam- 
ples of basic FK5 stars. The mean errors given in Ta- 
ble 3 refer to one 'mean' coordinate component. It is ob- 
tained as an rms average over £^ ia * and £^5, and over all 
the stars in the corresponding sample. The error budget 
for the 1535 basic FK5 stars is slightly fictious, since this 
sample contains double stars for which the FK6 provides 
in reality 'special' solutions instead of the 'direct' com- 
bination solutions discussed in this paper. Nevertheless, 
the results for this sample provide a valid indication for 
the overall accuracy of our combination method in the di- 
rect modes SI, LTP, and STP. The sample of 1202 basic 
FK5 stars contains 'apparently single objects' only. Most 
of these stars (878 objects) have direct solutions in the 
FK6. The error budget for these 878 basic FK5 stars in 
Part I of the FK6 is given in Wielen et al. (1999d). The 
error budgets for 3272 additional fundamental stars with 
direct solutions in the three modes are presented in Part 
III of the FK6 (Wielen et al. 2000a). 

Table 4 gives the error budget for the combination of 
the GC (Boss et al. 1937) with HIPPARCOS. The sample 
of GC stars are the 'full sample' of 29 717 GC stars ob- 
served by HIPPARCOS, and the 'subsample' of 11 737 GC 
stars with linear HIPPARCOS standard solutions. From 
Table 4 it is obvious that the original proper motions (iqc 
of the GC do not contribute significantly to the GC+HIP 
on average. However, for some brighter and well-observed 
stars, the accuracy of [iqc is much better than the rms 
value of e^.cc seen in Table 4. Furthermore, the proper 
motion /Zo(gc)Hj derived from the central positions of the 
GC and the HIPPARCOS Catalogue, has usually a quite 
reasonable accuracy for most of the GC stars. 

The typical gain in accuracy in the proper motions 
derived in the long-term prediction mode, relative to the 
original HIPPARCOS proper motions, is a factor of 4.6 for 
the basic fundamental stars in the FK6 = FK5+HIP, and 
a factor of 1.8 for the 11773 GC stars in the GC+HIP. 
This improvement in the LTP mode with respect to HIP- 
PARCOS is a consequence of the cosmic errors in the 
instantaneously measured HIPPARCOS proper motions. 
In contrast to the LTP mode, the short-term predictions 
(STP mode) do not differ so much from the HIPPARCOS 
solutions. However, in most cases we are more interested 
in the long-term averaged proper motion (LTP) or in the 
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Table 2. Results of the combination method for the star HIP 9884 = FK 74 = GC 2538 = a Ari for the various 
modes 







FK5 + 


HIPPARCOS 






GC + HIPPARCOS 




Quantity 




mean error S 


mean error 


a* 


mean error 


S 


mean error 


Input data: 


















A^f(Tf) or Ao;gc(Tgc) 


- 7.83 


± 12.52 


+ 96.95 


± 15.30 


- 2.64 


± 49.92 


+ 236.68 


± 32.88 


AflF or A/iGC 


+ 0.49 


-i- n /in 


i on 


± 0.38 


+ 1.69 


± 1.24 


3.98 


± 1.86 


Ui iQC 


1 Q47 84 




1929.73 




1 SQ9 (SO 




1890.30 




Aa;H(rH,i nd ) 


0.00 


± 0.77 


0.00 


± 0.54 


0.00 


± 0.77 


0.00 


± 0.54 


A/iH 


0.00 


± 1.01 


0.00 


± 0.77 


0.00 


± 1.01 


0.00 


± 0.77 


Tk.ind 


1991.26 




1991.51 




1991.26 




1QQ1 ^ 1 






49.48 


± 0.99 






49.48 


± 0.99 








17.28 




17.28 




17.28 




17.28 






2.91 




2.91 




2.91 




2.91 




A/Uo 


+ 0.18 


± 0.29 


- 1.57 


± 0.25 


+ 0.03 


± 0.51 


- 2.34 


± 0.32 


Results for the single-star mode SI: 


















Ax S i(T S i) 


- 0.18 


± 0.76 


+ 0.14 


± 0.54 


- 0.23 


± 0.77 


+ 0.08 


± 0.54 


A/jsi 


+ 0.23 


± 0.23 


- 1.34 


± 0.20 


+ 0.04 


± 0.42 


- 2.00 


± 0.29 


T S i 


1991.12 




1991.47 




iyy±.zD 




1991.51 




Apsi 


- 0.57 


± 0.95 






- 0.98 


± 0.95 






Results for the long-term prediction mode LTP: 














Axltp (Tltp ) 


- 5.14 


± 10.14 


+ 54.37 


± 11.46 


- 0.30 


± 16.34 


+ 51.26 


± 15.30 


A/JLTP 


+ 0.36 


± 0.31 


- 1.38 


± 0.27 


+ 0.28 


± 0.49 


- 2.37 


± 0.36 


Tltp 


1962.77 




1956.76 




1980.68 




1969.40 




Results for the short-term prediction mode STP: 














Axstp(Tstp) 


- 0.01 


± 0.77 


0.00 


± 0.54 


- 0.02 


± 0.77 


0.00 


± 0.54 


ApsTP 


0.00 


± 0.95 


- 0.08 


± 0.74 


- 0.03 


± 0.95 


- 0.15 


± 0.74 


Tstp 


1991.26 




1991.52 




1991.26 




1991.52 




ApsTP 


- 0.03 


± 0.99 






- 0.06 


± 0.99 







Units: mas, mas/year, or years 



single-mode result (SI), where the gain in accuracy is quite 
significant. 

9. Problems and applications 

The long-term predictions and the short-term predictions 
derived in the former sections are statistically the best as- 
trometric solution for a sample of 'apparently single-stars'. 
However, an unsatisfactory property of these statistically 
valid solutions is the fact that we are not able to make 
proper use of the available information on the individual 
behaviour of the stars. For example, we use the overall 
cosmic errors c M (p) and c x (p), no matter whether the ob- 
ject is a A/x binary or a single-star candidate (Wielen et 
al. 1999a). In principle, for each individual star, one would 
have to use 'conditioned correlation functions' which are 



based on the individually observed differences between the 
instantaneous measurements and the mean data (e.g. on 
/XFK5 — /xhip)- Unfortunately, the conditioned correlation 
functions, i.e. 'conditioned cosmic errors' in particular, are 
not available at present. 

We should also point out that presently the cosmic 
errors c x in position are much more uncertain than the 
cosmic errors c M in proper motion. While the typical values 
of are larger than the measuring errors e^fks, £^i.o, 
and e^Hip, the typical values of c x are nearly lost in the 
measuring errors of the ground-based data. 

Furthermore, our knowledge about the actual form of 
the correlation functions £(At), i](At), ((At) for A|=0 is 
presently still quite rudimentary. The simple example for 
the correlation functions presented in Sect. 3.6 of WPSA is 
mainly given for illustrating the general behaviour of these 
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functions. The application of this example to real data 
should be done very cautiously. Hence the real transition 
from the STP solution to the LTP solution (discussed in 
Sect. 5) is quantitatively not well-determined. 

Which of the three solutions offered (SI, LTP, STP) 
should be used in real applications ? For single-star candi- 
dates, the single-star mode should be adopted, although 
some of the single-star candidates may nevertheless be bi- 
naries. If the user is handling a sample of 'apparently sin- 
gle stars' (which usually contains single-star candidates, 
A/x binaries, and intermediate cases), then the LTP or 
STP solutions are recommended, depending on the corre- 
sponding epoch difference At = t — Th (with Th ~ 1991). 

The problems discussed above are especially severe for 
objects detected as Afi binaries (Wielen et al. 1999a). 
For Afi binaries the difference between the instantaneous 
proper motion (/imp) and the mean one (e.g., /xfks or 
fi ) is sometimes much larger than the cosmic error c M 
expected on average. In such a case the weight w^.ltp 
(Eq. (53)) of the HIPPARCOS proper motion (p mP = M2) 
is higher than appropriate for this individual object. The 
derived mean proper motion ^ltp is then biased towards 
the HIPPARCOS proper motion, and the derived mean 
error e^ltp is too small. In the case of extreme Afi bi- 
naries, it is better to adopt a properly weighted mean of 
/ii (e.g. /iFKs) and /io as a prediction for the mean proper 
motion fi m . The main problem in such a procedure is the 
unknown individual value of the cosmic error in the HIP- 
PARCOS position, which is higher than the value of c x 
expected on average. This value enters into the weight of 
/io, and hence into the predicted value of /i m . We shall 
discuss this problem in more detail in a subsequent paper. 
In any case, the LTP solutions (and, of course, the SI so- 
lutions) for A/i binaries are inherently the least accurate 
ones among the class of direct solutions, because of the 
disturbing double-star nature of these objects. 

10. Summary and outlook 

In this Paper II, we have derived and discussed an ap- 
propriate method to combine a ground-based astrometric 
catalogue (such as the FK5 or GC) with the HIPPAR- 
COS Catalogue, taking cosmic errors (due to undetected 
binaries) in the quasi-instantaneously measured data into 
account. The method leads to long-term predictions (LTP 
mode) and to short-term predictions (STP mode), which 
are the limiting cases of the general solution. The general 
solution is a smooth transition from the STP to the LTP 
mode. The case of single stars with no cosmic errors was 
already treated in Paper I (SI mode). In a subsequent pa- 
per, we shall present 'special solutions' for known double 
stars. 
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Table 3. Error budget for FK6 proper motions 



Typical mean errors of proper motions 
(in one component, averaged over //a* and fj,g; units: mas/year) 



Sample of stars: 


1535 FK 


1202 FK 




rms aver. 


median 


rms £iv6r. 


median 


SI: 


single- 


•star mode 






HIPPARCOS 










random 


0.82 


0.63 


0.68 


0.61 


FK5 










random 


0.76 


0.64 


0.77 


0.67 


system 


0.28 


0.25 


0.28 


0.25 


total 


0.81 


0.70 


0.83 


0.72 


Mo 










random 


0.53 


0.43 


0.54 


0.45 


system 


0.24 


0.23 


0.25 


0.23 


total 


0.58 


0.49 


0.59 


0.51 


FK6 = FK5+HIP 


0.35 


0.33 


0.35 


0.34 


i alio 01 1 1 1 1 i i \ i i v . v a j 


2.3 


1.9 


1.9 


1.8 












LTP: 


long-term prediction 




HIPPARCOS 










random 


0.82 


0.63 


0.68 


0.61 


cosmic (m m) 


2.15 


2.05 


2.13 


2.04 


total 


2.30 


2.14 


2.24 


2.13 


FK5 










random 


0.76 


0.64 


0.77 


0.67 


system 


0.28 


0.25 


0.28 


0.25 


total 


0.81 


0.70 


0.83 


0.72 


Mo 










random 


0.53 


0.43 


0.54 


0.45 


system 


0.24 


0.23 


0.25 


0.23 


cosmic (due to x) 


0.32 


0.30 


0.32 


0.30 


total 


0.66 


0.57 


0.68 


0.59 


FK6 = FK5+HIP 


0.49 


0.44 


0.49 


0.45 


ratio nf HTPPARPOS 


4.7 


4.9 


4.6 


4.7 


^ t~\ h l-v r\ V" V"f"\r"C 

LU r 1VU cllULa 










STP: 


short-term prediction 




HIPPARCOS 










random 


0.82 


0.63 


0.68 


0.61 


FK5 










random 


0.76 


0.64 


0.77 


0.67 


system 


0.28 


0.25 


0.28 


0.25 


cosmic (in n) 


2.15 


2.05 


2.13 


2.04 


total 


2.30 


2.16 


2.28 


2.16 


Mo 










random 


0.53 


0.43 


0.54 


0.45 


system 


0.24 


0.23 


0.25 


0.23 


cosmic (due to x) 


0.32 


0.30 


0.32 


0.30 


cosmic (in n) 


2.15 


2.05 


2.13 


2.04 


total 


2.25 


2.13 


2.23 


2.12 


FK6 = FK5+HIP 


0.68 


0.59 


0.62 


0.57 


ratio of HIPPARCOS 


1.2 


1.1 


1.1 


1.1 


to FK6 errors 
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Table 4. Error budget for GC+HIP proper motions 



Typical mean errors of proper motions 
(in one component, averaged over m«* and fig; units: mas/year) 



Sample of stars: 


29 717 GC 


11773 GC 




rms aver 


median rms aver. 


median 


SI: 


single- 


star mode 






HIPPARCOS 










random 


1.47 


0.73 


0.75 


0.69 


GC 










random 


10.57 


9.38 


8.59 


7.55 


system 


0.57 


0.49 


0.56 


0.48 


total 


10.59 


9.39 


8.61 


7.57 


Mo 










random 


1.78 


1.74 


1.42 


1.44 


system 


0.18 


0.12 


0.13 


0.12 


total 


1.79 


1.75 


1.43 


1.45 


GC+HIP 


0.72 


0.63 


0.62 


0.58 


lano Ol 1 1 1 1 1 J \ 1 \ \ A J .J 


2.0 


1.2 


1.2 


1.2 


to CJP4-HTP orrnrq 










LTP: 


long-term prediction 






HIPPARCOS 










random 


1.47 


0.73 


0.75 


0.69 


cosmic (in fx) 


1.75 


1.60 


1.80 


1.66 


total 


2.29 


1.76 


1.95 


1.80 


GC 










random 


10.57 


9.38 


8.59 


7.55 


system 


0.57 


0.49 


0.56 


0.48 


total 


10.59 


9.39 


8.61 


7.57 


Mo 










random 


1.78 


1.74 


1.42 


1.44 


system 


0.18 


0.12 


0.13 


0.12 


cosmic (due to x) 


0.12 


0.10 


0.12 


0.11 


total 


1.79 


1.75 


1.43 


1.45 


GC+HIP 


1.17 


1.15 


1.06 


1.05 


ratio Ol 1 1 1 1 1 J \ l\V AJ.J 


2.0 


1.5 


1.8 


1.7 


tr> PP-I-HTP prrnrs 
to v i v | 1 1 1 1 eiioib 










STP: 


short-term prediction 




HIPPARCOS 










random 


1.47 


0.73 


0.75 


0.69 


GC 










random 


10.57 


9.38 


8.59 


7.55 


system 


0.57 


0.49 


0.56 


0.48 


cosmic (in fi) 


1.75 


1.60 


1.80 


1.66 


total 


10.73 


9.53 


8.79 


7.75 


Mo 










random 


1.78 


1.74 


1.42 


1.44 


system 


0.18 


0.12 


0.13 


0.12 


cosmic (due to x) 


0.12 


0.10 


0.12 


0.11 


cosmic (in fj,) 


1.75 


1.60 


1.80 


1.66 


total 


2.51 


2.37 


2.30 


2.20 


GC+HIP 


0.82 


0.69 


0.69 


0.65 


ratio of HIPPARCOS 


1.8 


1.1 


1.1 


1.1 



to GC+HIP errors 



